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Long-Bone Allometry of Terrestrial Mammals and the 
Geometric-Shape and Elastic-Force Constraints of Bone Evolution 



A natural similarity in body dimensions of terrestrial animals noticed by 
ancient philosophers remains the main key to the problem of mammalian 
skeletal evolution with body mass explored in theoretical and experimental 
biology and tested by comparative zoologists. We discuss the long-standing 
problem of mammalian bone allometry commonly studied in terms of the 
so-called "geometric", "elastic", and "static stress" similarities by McMahon 
(1973, 1975a, 1975b). We revise the fundamental assumptions underlying 
these similarities and give new physical insights into geometric-shape and 
elastic-force constraints imposed on spatial evolution of mammalian long 
bones. 

A realistic description of an animal skeletal bone might be thought in 
terms of a hollow, irregular, curved structural rigid beam, whose linear di- 
mensions guarantee the avoidance of fracture caused by axial and non-axial 
external peak loads. Meanwhile, most mammalian skeletal bones manifest a 
geometric similarity that permits one to introduce a cylinder- shape approx- 
imation through the characteristic dimensions Dig and Lis related to bone 
diameter and length for a given bone i of a certain mammalian species s. 
This surmise can be verified experimentally with the help of bone allometry. 
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which provides data on mammalian allomctric exponents di and /j through 
the power-law scaling regressions Di^ = CigM'^^ and Lj^ = hisM^\ where body 
mass M M'i'°^'^^) is treated as an external mammalian parameter and 
and his are constants. The i-bone allometric exponents obey the cylindric- 
shape evolution constraint equation 2di -\- U — 1. This immediately follows 
from (i) the relation between hone mass Mis a-nd its volume Df^L^s {M^g — 
pDf^Lis , p is bone density), (ii) the aforegiven scaling regressions, and from 
(iii) McMahon's hypothesis that skeletal mass is linearly scaled to body mass. 
I.e., Mr ' oc Mis oc M. In what follows, we give an analysis of the bone 
dimension scahng in terms of the "overall-bone" (hereafter effective bone) 
mammalian allometric exponents d — n~^Ti"^^-^di and I — n^-'^Ef^^Zj, which 
obey the corresponding cylindric-shape constraint 2d-\-l — 1 resulting from 
the corresponding i-bone constraint averaged over n mammalian bones. Be- 
low we give statistical analysis for experimental data by Christiansen (1999b) 
on allometric exponents for mammalian long bones and show that it corrob- 
orates the cylindric-shape approximation very accurately. 

More sophisticated bone allometric studies provide information on evo- 
lution of real non-circular hollow bones with a cross section A^f^^ (cx M'^^, 
approximated by coaxial cylinders of length Lj^, with external Dis and inter- 
nal dis diameters and of area Ais = T:{D'jg — dlg)/A. As follows from Table 2 by 
Selker and Carter (1989) given for allometric exponents for Artiodactyl long 
bones, the circle-cross-section geometric similarity explicit in the constraint 

= 2di works well (with the accuracy amounting to 5%). Furthermore, 
similar to the internal bone diameter, the bone curvature, also shows a great 
deal of variation between different z-bones in the same animal and among dif- 
ferent s-animals, but does not show a significant influence from animal mass 
(Selker and Carter, 1989). We see that geometric similarity does not depend 
on local deviations of bone shape from the ideal cylinder form and can there- 
fore be castled into the cylindric-shape evolution constraint equation. It is 
worth noting that the geometric similarity by McMahon (1975b) corresponds 
to a particular case of the isometric solution of the cylindric-shape constraint 
given by dQ = Iq = 1/3 and oq = 2/3. 

Generally speaking, elastic-force mammalian similarity suggests that there 
exists a unique mechanism, attributed to the natural elasticity of the bones 
of adult animals, that provides an escape from the critical elastic deforma- 
tion (fracture) of bones under peak mechanical stresses. This similarity is 
ensured by a safety factor of about 2-4, independent of mammalian mass 
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(see Biewcncr. 1990 among others). Following Rashevsky (1948) and McMa- 
hon (1973), the elastic similarity can be introduced in explicit form on the 
basis of a mechanical analogy that takes place between a given bone (or an 
animal trunk) and a quasi-cylindrical, rigid beam. Elastic-force mammalian 
similarity, in a way, completes the geometric-shape similarity that eventu- 
ally provides predictions for the dimension-growth effective bone exponents 
d and / tested by long bone allometry. For the case of long bones, it seems 
plausible to approximate the long {Dis <C Lis) cylindrical beams by rigid 
rods for which the condition of elastic instability against axial loads was first 
established by Euler. This was given in terms of the critical elastic-buckling 
force F^ucki~ CTr^E^/js/Lf^ where E is the elastic modulus and I is — nDis^/QA 
is the second cross-area moment of inertia. A numerical factor c depends on 
the boundary conditions, i. e. on the way of application of the axial external 
forces at the rod ends. One can show (see e.g. Landau and Lifshits, 1989) 
that c — 1 and 4 for the cases of both the non-fixed (on the hinges) ends 
and fixed ends, respectively. The first case was cited by Hokkanen (1986) 
concerning with a clarification of the elastic similarity employed by McMa- 
hon and Rashevsky. Then, Carter and Spengler (1982, cited by Selker and 
Carter, 1989) reported that long bones in living animals are extremely rarely 
fractured by axial loads. This is due to the curved structure of bones, which 
in the case of in vivo axial loading introduce bending moments in their dia- 
physis, that create stresses of a greater amplitude than those caused by pure 
axial loading (Selker and Carter, 1989). This real case can also be related to 
Euler's solution with c = 1. McMahon (1973) extended these conditions for 
the case of the limbs of trees, when a free horizontal bending of the top from 
the trunk is caused by proper weight. This results in c = 1/4 and permits 
one to summarize all the buckling-deformation cases for the critical elastic 
force as 

7-)4 

elast ~ ^buckl ^ t2 ' \^) 
^is 

McMahon (1973) suggested additionally that all bone elastic forces for 
terrestrial mammals arc caused solely by animal weight {F^ucki^ 9 
is the gravitation constant) and therefore McMahon's elastic similarity for 
mammalian long bones can be introduced (see also discussion by Hokkanen, 
1986) through (i) the elastic-rod, critical-force and (ii) the cylindric-shape 
constraints, namely 



3 



a) M- 21^1 and b)2d + l^l. (2) 

Similarly to the aforegiven b constraint, the elastic-force a constraint in (2) 
follows from the allometric bone-dimension scaling applied to (1) further av- 
eraged over all long bones. The system of two equations provides the well 
known (McMahon, 1973) elastic model predictions for the bone diameter 
^{bucki) _ 3y'g^ bone length [^'"'^'''■^ = 1/4 exponents, along with the length- 
to-diameter [t^'^^^^ / do^'^^''^) dimension-bone exponent — 2/3, related 
to the allometric scaling Lig ^ D^^. The latter was in part corroborated 
by bone allometry by McMahon (1975b) who derived A^^^'* = 0.67, 0.52 
and 0.83 for families, respectively, Bovidae, Suids, and Cervids in the or- 
der Artiodactyla. A similar study by Alexander (1997) established d^^^^ 
= 0.34 and Zg^^^^ = 0.26, with aJ^^^^ = 0.76 for the family Bodivie. When 
a wider variety of phylogenetically and anatomically distinct samples were 
analyzed within a larger body size range, experimental allometric exponents 
exhibited rather the isometric scenario (Aq = 1) than the elastic model be- 
havior predicted by McMahon (1973, 1975a). This incompatibility with long- 
bone allometric observations was reported among others by Alexander et al. 
(1979a) and Biewener (1983). Furthermore, recent systematic investigations 
by Christiansen (1999a, 1999b) suggest that experimental observations of 
bone dimensions "support neither geometric nor elastic similarity, making 
both questionable as a means of explaining long-bone scaling in terrestrial 
mammals". Moreover, according to criticism by Economos (1983) shared by 
Christiansen (1999b), no satisfactory explanation for any power-law scaling 
observed in mammalian allometry can be expected. Prom the physical point 
of view, however, there is no doubt that long bones possess simultaneously 
cylindrical-type-shape and elastic-critical-force similarities. We therefore re- 
vise the underlying hypotheses of the elastic model by McMahon (1973) to 
describe the elastic-force {F^Z^'^ oc M and F^^} = ^et^) and geometric- 
shape [Ms^^^^^ ^ M) similarities in terms of constraints given in (2). 

First, under locomotion peak dynamic conditions, skeletal clastic forces 
of animals may exceed those driven by static body weight (Hokkancn, 1986), 
and therefore the peak stresses in bones are due to muscle contractions, rather 
than to gravity (Carter et al., 1980, Biewener, 1982, Rubin and Lanyon, 1984, 
Selker and Carter, 1989, and Biewener, 1991). This implies that the peak 
elastic forces in (1) should be substituted by those of muscle subsystem (mus- 
cles, tendons and ligaments), which stores and return elastic energy (Farley 
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et al. 1993), i.e., F^'^^^ = F^^\ Furthermore, after Rubin (1984) many 
studies provide strong evidence that the maximum muscle strains recorded 
in rigorous activities of animals are independent of body mass, and thus 
muscle-induced bone stresses are Fg^st^ /^m M^, where Aj^ is the cross- 
section area of muscle fibers. This, in turn, exhibits scaling to body mass 
through the muscle-area allometric exponent defined through the scal- 
ing law A„i oc M""^. As a consequence, the elastic-buckling-force constraint 
should be modified and changed in (2) for 4d — 21 = a, with a = a^- 

Second, in view of the observation by Prange et al. (1979) of nonlinear 
evolution of skeletal mass with body mass M, Hokkanen (1986) noted that 
the corresponding McMahon's hypothesis (m]*'^^'^ ^ M) makes the b con- 
straint inaccurate in (2). With taking into account the bone-mass-to-body- 
mass scaling law Mis = pDfgLis oc one obtains the following modified 
cylindrical-shape bone constraint for terrestrial mammals: 2di-\-li = h^. Mod- 
ified in such a way McMahon's constraint results in a) 4d — 21 — a and b) 
2d-\- 1 — b that provides a new prediction due to elastic-buckling similarity, 
namely 

^ibucki) ^ a(-P) + 26(-P) ^ 2&(-P) - a(-P) 

8 4 ■ ^ ^ 

This prediction is tested below on the basis of the experimental data on 
muscle 

(exp) g^j^^ bone-mass h^^^^^ mammalian effective exponents. 
Third, after Alexander et al. (1979b) it has been widely recognized that 
the external critical loads acting parallel to the axis of the bone shaft (di- 
aphysis) are often fewer that those applied in the perpendicular direction, 
although the corresponding peak stresses are of the same order of magnitude 
(see Rubin and Lanyon, 1984). Within the context of theory of elasticity, 
this experimental finding can be treated in terms of the thermodynamic in- 
stability of an ?'-s-bone approximated by a convex shell of characteristic size 
Lis with fixed ends subjected to deformation of flexure Hig. The free energy 
of the shell (a difference between the elastic energy and the work of defor- 
mation produced by some uniform external peak pressure p('=''**)) shows its 
instability when j/^'**) E{h\^^^ / LisY (Pogorelov, 1960). Adopting for a 
critical condition of bone fracture H^T^^ ^ Dig one has 



Tp(crit) rp(crit) 
p{crit) ^ buckl ^ bend ^ ^ is 
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A scaling of the critical elastic-buckling force to bone dimensions, that follows 
from (4), is given in (1), and that for the elastic-bending case is 



j-,(max) _ rn(crU) rpicrit) j^t^ ^ ^(crit) r j n (^\ 

^elast - ^bend ^ ^tors ^ ^-j— > ^ buckl ^Or Li, > Vi,. (5j 

Besides the critical force -F^m? ^ perpendicular load that a rigid beam 

can withstand without breaking that was already cited by Hokkanen (1986) 
and by Selker and Carter (1989, see (11) and (4), respectively), we have 
included a corresponding torsional force F^^"g^\ An estimate for the latter 
immediately follows from the critical torsional angle of a slightly bent and 
twirled rod considered by Landau and Lifshits (1989). For the critical tor- 
sional moment along the rod axis, one therefore has F^^f^'Dig — 8.98EIis/Lis 
that results in the scaling given in (5). 

Although the importance of bending and torsional critical loads in the 
production of peak bone stresses is well established (Rubin and Lanyon, 1982, 
Biewener, 1982, Biewener and Taylor, 1986, Selker and Carter 1989), we 
have demonstrated a natural property of elasticity of bones, similar to the 
geometric-shape similarity, in a certain way does not depend on details of 
the local bone geometry and bone-end boundary conditions, i.e. all possible 
elastic-force scaling (similarities) follow from (4). The elastic-force bending- 
torsional scaling given in (5) provides the corresponding critical-force con- 
straint 3d — I — a, that in combination with the modified cylindrical-shape 
constraint 2d + I — b results in the relevant bending-torsional criterium for 
effective-bone evolution, namely 

^ibend) ^ a +0 ^ 60 Za 

5 5 

As seen from (6), within the context of the hypotheses adopted by McMa- 
hon (1973) one has 4^'"'^^ = 2/5, ^''""'^^ = 1/5, and A^"''^ = 1/2 that was 
denominated as a static stress similarity (McMahon (1975b). Remarkable 
that both the elastic-force criteria, given in (3) and (6) show a consistence 
with the isometric solution. For this specific case, the muscle subsystem de- 
velops independently (69 = 1) and isometrically {oq = 2/3) and both the 
critical elastic forces scale to mass as i^jj"^^^-* a M^/'^ regardless of the under- 
lying bone-structure protecting mechanism against critical buckling, bending 
or torsion deformations (see similar discussion by Selker and Carter, 1989). 
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Let us test the predictions for the bone-diameter Sp'^^'^') and /(p^^'^) growth 
exponents given in (3) and (6) with those known from the one-scale (overall 
small and large mammals) bone allometry and those available in the liter- 
ature data on a'^'^^P-* and ^('^^p). For the case of a constraint, an estimate 
for the overall mammalian data, including birds, on muscle-area allomet- 
ric exponent obtained among others by Alexander (1977), Alexander et al, 
(1981), and Pollock and Shadwick (1994), i.e., a'm^^ = 0.77 - 0.83, was pro- 
posed by Garcia (2001). The most recent data on the i-bone-mass allometric 
mammalian exponents ^•'^^^•'by Christiansen (2002) are presented in the last 
column of Table 1, and the relevant long-bone mammalian exponent ^(^^p) 
(= 'Ej^ibf^'^^ /A) can be introduced with accounting for statistical error due 
to distinct regressions, approximately, by 6('^^p) — 1.0 — 1.1. In Fig.l we give 
a comparative analysis for the distinct bone-evolution scenarios predicted by 
the elastic-force buckling (area 1) and bending-torsional criteria (area 2), on 
the basis of (3) and (6), respectively, with the help of the aforegiven experi- 
mental data a^^'^P^ = aln"^'^ and b^^^^\ The modified and improved cylindric- 
shape and elastic-force constraints are shown by solid lines for the limiting 
data a*^^^P^ and ^('^^p) and the dashed lines and open points correspond to the 
simplified constraints by McMahon. In general, one can see that the long- 
bone allometry data on mammalian effective (all-bone-averaged) exponents 
by Biewener (1983), Bertram and Biewener (1992, reestimated by Garcia, 
2001), and Ghristiansen (1999a) are close to both the elastic-force criteria. 
An exception although should be made for the pioneering data by Alexander 
et al. (1979a), where dispersion due to phylogenetic spectrum of terrestrial 
mammals was not reduced (see discussion by Ghristiansen, 1999b). 

Besides the effective-bone mammahan evolution observed from Fig.l, we 
give a statistical analysis of the similarity constraints in Table 1 with the 
help of the one-scale^ i-bone allometric data on the dimension exponents 
df^'^\ if^"^"^ and the i-bone-mass exponents bf^"^^ derived by Ghristiansen 
(1999a, 1999b and 2002, respectively) within four decades of mammalian 
mass {1-lO^kg) through the two different one-scale regressions. As seen from 
Table 1, the cylindrical-shape approximation for the geometric-shape simi- 
larity through the constraint 2(1^^"^^+ if^"^^ = bf^"^^ is experimentally justified 
with a good level of accuracy by the data on bf^^^ derived by both the meth- 
ods for all n = 1,2. ..6 long bones, except for the case of the ulna. A high 

^Extended statistical analysis of the long-bonc allometric data scaled for small and 
large mammals by Christiansen (1999a, 1999b) was given by Kokshenev (2003). 
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precision is established for the bone-mass allometric exponent b^*^^"^^ = 1.0-1.1 
predicted on the basis of the cylindrical- shape constraint b^^^'''^ — 2 o?(®^p)-|- 
/(cxp) -^itiiin 4-bone averaging, with b'^P'''^'^') = 1.03-1.06. 

The elastic-force similarities are tested through the predicted effective- 
bone mammalian exponents a^'^^^\ related to the muscle-area exponents, 
found as certain predictions given by the elastic-force buckling and bending- 
torsional constraints , respectively: a(^"'=^0 = 4^*^'^''^) — 2/'^'^''^) and a(*'="'^) = 
2^(cxp) _^(cxp)^ follows from the analysis elaborated in Table 1, the overall- 
bone predictions for a^^'^'^^^^ and a^^^'^'^) are, respectively, wedged between 0.89- 
0.97 and 0.81-0.87. With taking into account that is hmited by 0.83, 

we infer that the predicted data based on (i) that the elastic buckling force 
given in (1) and those for (ii) long bones fibula and ulna (treated by the RMA 
regression method) contradict to muscle-area allometry and therefore should 
be excluded from the critical-force evolution scenario. The first statement 
is discussed in (4) and the second is in line with the notes by Christiansen 
(1999a, 1999b) that " too thin u/no" and greatly reduced fibula " do not have 
"any importance in support of body mass" . In the special case of carnivore's 
ulna, however, provides a substantial forelimb support (Christiansen, 1999a). 

As follows from analysis given in Table 1, the overall-bone {LSR method) 
prediction for the muscle-area exponent is a^^^"'') = 0.81-0.83, which coin- 
cides with the critical-muscle-force exponent a^^^ = 0.81-0.83 provided by 
muscle-fiber allometry data by Pollock and Shadwick (1994) reanalyzed re- 
cently by Kokshenev (2003). This finding implies that the statistical analysis 
of the elastic-force criteria based on the equation -Fg;™^f "* = is related 

to the peak-area muscles and should therefore exclude muscle fibers with 
relatively small cross-sections. This can be exemplified by common digi- 
tal extensors which have the lowest "maximum" area and exhibit isotropic 
evolution with a™^^'' ~ = 2/3, as shown in Fig.l by Kokshenev (2003). 
Finally, let us estimate the aforementioned length-to-diameter effective bone 
exponent A^'"^'''^) = {'ib^''^'^'^ - 2a^''^'^^) / {b'^"^'^^ + a^P^^'^^) predicted with 4-bone 
exponents b^^'^'^ = 1.03-1.06 and a^^'^'^ = 0.81-0.83 obtained in Table 1. 
This provides A'-*'^"'^-' = 0.77-0.83 that corroborates the allometry exponent 
^(cxp) ^ 0.781-0.810 derived by Kokshenev (2003) from the relevant data by 
Christiansen (1999b). 

To summarize, we have revised McMahon's elastic similarities that have 
long been a controversial subject of intensive study in the last three decades, 
especially in the long-bone allometry for terrestrial mammals. The elastic- 
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force and cylindrical-shape bone evolution constraint equations were not sup- 
ported experimentally, and we have therefore reconsidered the basic hypothe- 
ses of the approach by McMahon (1973). As a matter of fact, exploration of 
the long bone allometric data was commonly given within the framework of 
model simplifications such as (i) the skeletal subsystem of animals operates 
separately from the muscle subsystem and (ii) the bone stresses are induced 
solely by gravity (Rashevsky, 1948, McMahon, 1973). These contradict to 
the observed evolution of bone mass and ignores a role of muscle contrac- 
tions in formation of peak skeletal stresses. McMahon's description of the 
cylindrical-shape and elastic-force similarities implicit in the bone dimen- 
sional constraints is revised, modified and extended in view of the up-to-date 
knowledge on the bone-mass and muscle-area scalings and on domination of 
the bending-induced deformations. 

Unlike the case of McMahon's (1973) elastic similarity model, the elastic- 
force (buckling-deformation) and cylindrical-shape similarities modified by 
muscle-area and bone-mass evolution (shown by shaded area 1 in Fig.l) be- 
come closer to almost all cited experimental data. At the same time, the 
most accurate data by Christiansen (1999a, 1999b) conflicts with the mod- 
ified elastic-buckling bone evolution scenario. Further improvement of the 
elastic-force similarity by critical bending and torsional deformations made 
it eventually observable (see shaded area 2 in Fig.l). Statistical analysis 
of the allometric data by Christiansen (1999b, 2002) corroborates, at least 
in the case of the standard LSR method, the proposed description of the 
geometric-shape and elastic-force mammalian similarities given in terms of 
the corresponding constraints within the cylindrical-shape and elastic-rod 
model approximations. 

The authors are indebted to R. McNeill Alexander for sending the detailed 
experimental data. Thanks are due to J. Tyson for valuable criticism, which 
improved substantially a presentation of the manuscript. Financial support 
by the CNPq (V.B.K.) is also acknowledged. 
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/(cxp) 
''i 




a-i 


k 


r,(cxp) 








.3860 


.3109 


0.92 


0.85 


1.08 


1.083 








.3599 


.3089 


0.82 


0.77 


1.03 


1.071 








.3654 


.2767 


0.91 


0.82 


1.01 


0.998 








.4014 


.3210 


0.96 


0.88 


1.12 


1.101 








.3681 


.2430 


0.99 


0.86 


0.98 










.4579 


.3177 


1.20 


1.06 


1.23 


1.101 








^(cxp) 


/(cxp) 


a 


a 


b 


5(cxp) 








.3782 


.3044 


0.90 


0.83 


1.06 


1.063 








.3762 


.2921 


0.92 


0.84 


1.04 










.3898 


.2964 


0.97 


0.87 


1.08 







Table 1. Analysis of the bone evolution mechanism through observation 

of the mammalian long-bone similarity on the basis of the elastic-force and 
cylindrical-shape evolution constraint equations. Experimental data on the 
mammalian bone diameter df'^'^\ length if'^'^^ allometric exponents for 6 
long bones obtained by the least square regression (LSR) and the reduced 
major axis (RMA) methods are taken from Table 2 by Christiansen (1999a). 
The RMA data for the bone- mass allometry exponent bf^'^\= b^^^'^^^) is 
taken from Table 2 by Christiansen (2002). The LSR data is estimated here, 
approximately, through the relation b* = where r^^^^^^ is the 
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corresponding correlation coefficient by Christiansen (2002). The muscle- 
and the bone-mass hi exponents are predicted with the help of the 
elastic-force ^ucAi/m^f-deformation (oj = 4(ij — 2/j) and ftent/m^'-deformation 
{tti = 3di — li) constraint equations, and cylindrical-shape constraint equation 
{bi — 2di+ li). The effective n-bone mammalian allometric exponents are 
given through the mean values of n corresponding exponents, e.g., a — 
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FIGURE CAPTURES 



Fig.l. Mammalian bone-dimension evolution diagram: diameter against 
length. Solid lines arc given by the elastic-force a-constraint buckling {4d — 
21 = a) and bending {3d — I = a) equations, and for the cylindrical-shape 
constraint {2d + I = b) equation. Dashed lines - the same for the case of 
a = b = 1. Points: {open squares) 0, 1 and 2 correspond, respectively, to the 
isometric (1/3, 1/3), the elastic similarity (3/8, 1/4), and to the stress simi- 
larity (2/5, 1/5) models by McMahon (1975b); closed circles symbolized by 
A'79, B'83, B'92 and C'99 arc due to one-scale, overall-bone- averaged, least- 
squarc-rcgression data on the mammalian long-bone allometry by Alexander 
et al. (1979a), Biewener (1983), Bertram & Biewener (1992, reestimated 
by Garcia, 2001), and Christiansen (1999a), respectively. Statistical error is 
shown by bars. Dashes areas 1 and 2 correspond to the proposed criteria 
given by the predictions d^^'^^^\ ^(''"'^'=0 and d^^^"'^\ litend) respectively, (3) 

and (6) for the case of the overall muscle areas with am^^"* = 0.77 — 0.83 (by 
Garcia, 2001) and bone masses with 6(^^P) = 1.0-1.1 (from Table 1). 
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